Scaling in the Fan of an Unconventional Quantum Critical Point 
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We present results of extensive finite-temperature Quantum Monte Carlo simulations on a SU(2) 
symmetric 5=1/2 quantum antiferromagnet with a four-spin interaction [Sandvik, Phys. Rev. 
Lett. 98, 227202 (2007)]. Our simulations, which are free of the sign-problem and carried out 
on lattices containing in excess of 1.6 x 10 4 spins, indicate that the four-spin interaction destroys 
the Neel order at an unconventional z = 1 quantum critical point, producing a valence-bond solid 
paramagnet. Our results are consistent with the 'deconfined quantum criticality' scenario. 



Research into the possible ground states of SU(2) sym- 
metric quantum antiferromagnets has thrived over the 
last two decades, motivated to a large extent by the un- 
doped parent compounds of the cuprate superconductors. 
In these materials, the Cu sites can be well described as 
S = 1/2 spins on a two-dimensional (2D) square lattice 
that interact with an anti-ferromagnetic exchange, the 
archetypal model for which is the Heisenberg model. By 
now, it is well established [l| that the ground state of this 
model with nearest-neighbor interaction has Neel order 
that spontaneously breaks the SU(2) symmetry. Two 
logical questions immediately arise: What possible para- 
magnetic ground states can be reached by tuning compet- 
ing interactions that destroy the Neel state? Are there 
universal quantum-critical points (QCP) that separate 
these paramagnets from the Neel phase? 

An answer to the first question is to disorder the Neel 
state by the proliferation of topological defects in the 
Neel order parameter It was shown by Read and 
Sachdev 0] that the condensation of these defects in 
the presence of quantum Berry phases results in a four- 
fold degenerate paramagnetic ground state, which breaks 
square-lattice symmetry due to the formation of a crystal 
of valence bonds - a valence-bond solid (VBS) phase. An 
answer to the second question was posed in recent work 
by Senthil et al. [4| , where the possibility of a direct con- 
tinuous Neel-to-VBS transition was proposed. The natu- 
ral field theoretic description of this 'deconfined quantum 
critical point' is written in terms of certain fractionalized 
fields that are confined on either side of the QCP and 
become 'deconfined' precisely at the critical point. As 
is familiar from the general study of QCPs, these frac- 
tional excitations are expected to influence the physics in 
a large fan-shaped region that extends above the critical 
point at finite-T Q (see Fig. [J). 

It is clearly of great interest to find models that harbor 
a direct Neel- VBS QCP and that can be studied with- 
out approximation on large lattices. Currently, the best 
candidate is the 'JQ' model, introduced by Sandvik 
which is an S = 1/2, SU(2) invariant antiferromagnet 
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FIG. 1: (color online) Schematic of the proposed T-J/Q 
phase diagram of the JQ model. Large-scale finite-T sim- 
ulations presented here substantiate the following: (i) The 
T = Neel order present for J/Q S> 1 is destroyed at a QCP 
(J/Q » 0.038); (ii) In the 'quantum critical fan', there is scal- 
ing behavior characteristic of a z = 1 QCP; (iii) An accurate 
estimate of the scaling dimension of the Neel field establishes 
that this transition is not in the 0(3) universality class; and, 
(iv) The paramagnetic ground state for sufficiently small J/Q 
is a VBS. In the QMC basis used here, the region with Q < 
is sign problematic. Fhe inset shows how the frustrating Q 
term is written in terms of bonds on a plaquette. 



with a frustrating four-spin interaction, 
H JQ = J^SrS.-Q ^(Si-Sj-ijCSfc-Sj-i), (1) 



(ijki) 



where indices are arranged as in the inset of Fig. [T] Using 
a T = projector Quantum Monte Carlo (QMC) method 
on lattices sizes up to 32 x 32 [f| , Sandvik showed that the 
four-spin interaction destroys Neel order and produces a 
VBS phase at J/Q ~ 0.04. Close to this critical value of 
J/Q, scaling in the spin and dimer correlation functions 
suggests a continuous transition with anomalous dimen- 
sions of the Neel and VBS order parameters equal, with 
a common value r/ = 0.26(3). In this Letter, we explore 
the candidate Neel- VBS QCP in the full T - J/Q phase 
diagram on large lattices using a complementary finite-T 
QMC technique, the Stochastic Series Expansion (SSE) 
method with directed loops 7]. The SSE QMC allows 
access to the physically important quantum critical fan 
(see Fig. [J), and admits high-accuracy estimates for the 
spin stiffness, ps, and the uniform susceptibility, Xu- The 
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FIG. 2: (color online) T — > converged Neel (main) and 
VBS (inset) order parameters as a function of 1/L. Dashed 
lines are quadratic fits that illustrate the finite condensate in 
the ordered phases. The solid (red) line is a fit to the form 
y — c\x C2 (illustrated for J/Q = 0.040), where C2 = z + nn 
is expected at the critical coupling. In fitting to the nine L 
values for each J/Q, we find a minimum in the chi-squared 
value (per degree of freedom) of 3.1 for J/Q — 0.040, with 
C2 ~ 1.35(1). For J/Q = 0.038, the chi-squared value is 3.9, 
with C2 ~ 1.37(1). All other J/Q produce much larger chi- 
square (greater than 10). 



scaling of these observables provides strong evidence for 
a continuous z = 1 transition in the JQ model. 

Basis and Sign of Matrix Elements: A priori, it is un- 
clear that SSE simulations of Hjq are free of the notori- 
ous sign-problem: a fluctuating sign in the weights used 
in the QMC sampling. In the SSE, finding an orthog- 
onal basis in which all off-diagonal matrix elements of 
the Hamiltonian are non-positive solves the sign-problem. 
A simple unitary transformation on the S z basis (a tt- 
rotation about the z-axis on one sub-lattice) results in a 
new basis in which, for J, Q > 0, all off-diagonal matrix 
elements of Hjq are non-positive, allowing sign-problem 
free simulations (Fig.[T|). We note that this non-positivity 
condition is also the main ingredient in the proof of 
the Marshall sign theorem, allowing us to infer that the 
ground state of -Hjq for J, Q > must be a spin-singlet. 
As shown below, this singlet state changes from Neel at 
Q < J to VBS at J < Q. 

Numerical Results: Using the SSE QMC, we studied 
various physical observables in the JQ model on finite- 
size lattices of linear dimension L (with number of spins 
A'spin = L 2 ). Particular attention was paid to the scal- 
ing of the spin stiffness p s = d 2 E /d<p 2 (E is the energy 
and <f> is a twist in the boundary conditions) and the uni- 
form spin susceptibility Xu = ((52i Sf) 2 } /TN sp \ n . In the 
S z basis used here, it is easy to measure the correlation 
functions C§(r,r) = (S z (r, t)S z (0, 0)) and C$(r,r) = 
([S*(r,T)S*(r + x,T)][S*(0,0)S*(x,0)]>. While C§ is the 
correlation function of the Neel order parameter, the VBS 
order is indicated by Cy, which is the correlation function 
of the composite operator S z (r)S z (r + x), receiving con- 
tribution from both the standard VBS order parameter 
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FIG. 3: (color online) Criticality of the Neel field at J = 0.038: 
collapse of the Neel structure factor (Sjv) and susceptibility 
(xn) with z — 1 and T]n = 0.35, determining the universal 
functions Xs(a;) and X x (x) (up to non-universal scale factors 
on the x and y axes). The only fit parameter for both Sn and 
Xjv is vn, the anomalous dimension of the Neel field. 



S(r) ■ S(r + x) as well as the traceless symmetric ten- 
sor constructed from S l (r)S 3 (r + x). Structure factors 
for the Neel and VBS phases are constructed from these 
correlation functions by Fourier transformation at equal 

time, 5 NjV [q] = ^ r [ ex P( _ic l ' r ) C< N,v( r ' T = °)]/^spm, 
from which the order parameters are defined at the ob- 
served ordering wavevectors: m 2 ^ v = S l N,v[qN,v]/A r sp i n . 
Zero-frequency susceptibilities (xn and Xv) ar e con- 
structed by integrating over all r and Fourier transform- 
ing in space to the ordering vectors, qN,v- 

Examination of the full q-dependent structure factors 
indicate the presence of sharp ordering wavevectors in 
<Sn[<1n = (tt, 7r)] for large J/Q and -Svfqv = {k, 0) or 
(0, 7r)] (the latter in the case where the correlator is mea- 



sured with y) for large Q/J 
VBS phases observed in Ref. 



] , confirming the Neel and 
|. As shown in Fig.H T -> 



converged data scales convincingly to a non-zero value 
for at J/Q = 0.1 and for m\ at J/Q = 0.01. The 
critical coupling appears to occur between J c « 0.038 
and 0.040 (we set Q = 1 fixed throughout), such that 
as J c is approached from above (below) the extrapo- 
lated Neel (VBS) order parameter is suppressed. Very 
near J c , both order parameters vanish within our er- 
ror bars, while a power law with no y-intercept fits the 
Neel data with high accuracy. More specifically, at J c , 
scaling arguments require Sn oc L 1 ~ r,N Xs(L z T/c) and 
Xn oc L 2 ~ VN H x (L z T/c), with t]n the anomalous dimen- 
sion of the Neel field. In Fig. [3l we verify this scaling 
behavior and determine the universal functions Xs and 
X x . Both analyses illustrated in Figs. (I2|.'il) give a consis- 
tent estimate of r]N ~ 0.35(3). This value is larger than 
the result of tjn ~ 0.26(3) from Ref. [fj. While the exact 
source of this discrepancy is unclear due to the entirely 
different methods used to extract the exponents, we note 
that (a) our analysis does not involve extra fit param- 
eters from the inclusion of sub-leading corrections, and 
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FIG. 4: (color online) Finite- T uniform susceptibility, for a 
L — 128 system near J c . Error bars are much smaller than 
the symbol size. For the region 0.08 < T < 0.18, the data is 
highly linear, and a straight-line fit for J/Q — 0.038 (shown) 
intercepts the origin within error bars. Intercepts of straight- 
line fits for all data sets are in the inset. From the slope of 
the linear-T behavior we obtain A x /c 2 = 0.0412(2). 
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FIG. 5: (color online) Zoom-in of the spin stiffness and sus- 
ceptibility close to the expected critical point, taken for sim- 
ulations cells of size L = 1/T. Data in the inset is scaled to 
get the best collapse for the largest system sizes, which occurs 
for J c ~ 0.039(1) and v ~ 0.68(4). 



(b) the collapse of both Sn and xn takes place over two 
and a half orders of magnitude of LT with only one com- 
mon fit parameter, tjn; both facts give us confidence in 
our estimate. The critical scaling of C v is more compli- 
cated; due to the aforementioned mixing-in of two order 
parameters, Cy is expected to receive two independent 
power-law contributions. Indeed, it is difficult to disen- 
tangle these individual contributions on the limited range 
of lattices sizes available, precluding us from verifying the 
proposal [|| that rjN — rjy. 

We now turn to an analysis of the scaling properties 
of Xu and p s in the hypothesized quantum critical fan 
region of Fig[T] Xu and p s , being susceptibilities of con- 
served quantities have no anomalous scaling dimension, 
and hence at finite-T and L in the proximity of a scale- 
invariant critical point, assuming hyper-scaling: 

Ps(T,L,J) = ^Y^f, gL y^, (2) 

Xu(T,L,J) = ^z(¥f,gLV»y (3) 

where g cx ( J — J c )/J c - At criticality (g = 0), it is easy 
to see that Y(x — ► 0, 0) = A p /x and Z(x — > oo,0) = 
A x x d / Z , where Y(x,y) and Z(x, y) are universal scal- 
ing functions and A x , A p are universal amplitudes of the 
quantum critical point; c is a non-universal velocity. 

At criticality and L — > oo, one can show from Eq. ([3]) 
that Xu — -rhT d l z ~ x \ i.e. for a z = 1 transition, Xu 
should be T-linear and have a zero intercept on the y- 
axis at T = 10]. In Fig. 0J X u data for an L = 128 
system is presented. Within our error bars, this data is 
L — > oo converged for the region of T shown; at smaller 
T the finite-size gap causes an exponential reduction in 



Xu- The inset shows how the extracted value of the y- 
intercept, a (from a fit to the form a+bT), changes sign as 
the coupling is tuned, consistent with 0.036 < J c < 0.040 
and demonstrating to high precision the z — 1 scaling. 

Turning to study Eqs. further, one may hold the 
first argument of the universal functions fixed by setting 
L = 1/T (assuming z — 1 as indicated above). In order 
to achieve this, we performed extensive simulations on 
lattices sizes up to L = 1/T = 64, illustrated in Fig. [5] 
According to Eqs. ()2I3[1 . data curves for Lp s and Lx u 
plotted versus J should show a crossing point with dif- 
ferent L precisely at J c . We find that for relatively large 
sizes (32 < L < 64) the crossing point converges quickly 
in the interval 0.038 < J < 0.040. The insets show the 
data collapse when the x-axis is re-scaled to gL x l v (with 
v = 0.68). We note that with the inclusion of small 
sub- leading corrections (of the form a u /L u ), the cross- 
ing point and data collapse of p s and Xu can be made 
consistent, at the expense of two more fit parameters, 
even for much smaller system sizes than illustrated 0. 
In contrast to the U(l) symmetric JK model [Hj], where 
the absence of a T-linear Xu and a crossing in the data for 
p s L cast doubt on its interpretation as a z — 1 QCP, the 
present data for this SU(2) symmetric model gives strong 
support for a z = 1 QCP between 0.038 < J < 0.040. 

Finally, we hold the second argument of the scaling 
functions [Eqs. (|2I3[) ] constant by tuning the system to 
g = 0. One then expects a data collapse for p s /T and 
Lxu when they are plotted as a function of L Z T (with z = 
1). Fig. [6] shows this collapse for simulations carried out 
with extremely anisotropic arguments LT, varying over 
almost three orders of magnitude. There is an excellent 
data collapse over 8 orders of magnitude of the range of 
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FIG. 6: Scaling of \u and p s at J = 0.038 « J c , with z = 1 
and d = 2. These plots are the universal functions Y(x, 0) and 
Z(x, 0) up to the non-universal scale factor c on the i-axis. 
The expected asymptotes (see text) are plotted as dashed 
lines Y(x -> 0, 0) = A P /x and Z(sc -» oo, 0) = A x x d/z . From 
fits to the data, we find A x /c 2 = 0.041(4) and A p c = 0.37(3), 
allowing us to estimate a universal model-independent num- 
ber associated with the QCP, A P ^[A^ ~ 0.075(4). 



the universal functions, with no fit parameters. This data 
together with that in Fig. |4] provide our most striking 
evidence for the existence of a QCP with z = 1 in the 
proximity of J/Q as 0.038. 

Discussion: In this paper we have presented extensive 
data for the SU(2) symmetric JQ model which indicates 
that the Neel order (present when J ^> Q) is destroyed 
at a continuous quantum transition as Q is increased . 
In the finite- T quantum critical fan above this QCP, scal- 
ing behavior is found that confirms the dynamic scaling 
exponent z = 1 to high accuracy. The anomalous dimen- 
sion of the Neel field at this transition is determined to be 
?7n ~ 0.35(3), almost an order of magnitude more than 
its value of 0.038 12] for a conventional 0(3) transition. 
For sufficiently large values of Q we find that the system 
enters a spin-gapped phase with VBS order. To the ac- 
curacy of our simulations, our results are fully consistent 
with a direct continuous QCP between the Neel and VBS 
phases, with a critical coupling between J/Q as 0.038 
and J/Q w 0.040. Although our finite size study can- 
not categorically rule out a weak first-order transition, 
we have found no evidence for double-peaked distribu- 
tions, indicating an absence of this sort of first-order be- 
havior on the relatively large length scales studied here. 
It is interesting to compare our results to the only the- 
ory currently available for a continuous Neel- VBS tran- 
sition: the deconfined quantum criticality scenario 0] , in 



which the Neel- VBS transition is described by the non- 
compact CP 1 field theory. All of the qualitative observa- 
tions above, including an unusually large 7?n [H| agree 
with the predictions of this theory. Indeed, our esti- 
mate of 7?n ~ 0.35 [Fig. [3] is in remarkable numerical 
agreement with a recent field-theoretic computation 14 1 
of this quantity, which finds ??n = 0.3381. With regard 
to other detailed quantitative comparisons, we have pro- 
vided the first step by computing many universal quanti- 
ties, X x (x), X s (x), Y(ar,0), Z(x,0) and A PV ^ « 0.075 
[Fig. [5] in the JQ model. Analogous computations in the 
CP 1 model, although currently unavailable [lj| are highly 
desirable to further demonstrate that the JQ model re- 
alizes this new and exotic class of quantum criticality. 
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